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Abstract 

A theoretical study on the phase behaviors of a polymer solution confined between two coaxial 
cylindrical walls is presented. For the case of a neutral inner cylinder, the spinodal point derived 
through the Gaussian fluctuation theory is confinement-independent because of the existence of 
a free dimension in the system. The kinetic analysis indicates that the fluctuation modes always 
have a component of a plane wave along the axial direction, which can lead to the formation of 
a periodic-like concentration pattern. On the other hand, the equilibrium structure of the system 
is obtained by using the self-consistent mean-field theory (SCMFT) and the interplay between 
the "wetting" phenomenon and the phase separation is observed by modifying the property of the 
inner cylinderical wall. 
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I. INTRODUCTION 



Rigorously speaking, the classical thermodynamics should be always used to deal with 
the systems of infinite particles in an inifinite space. For a real finite system with a large 
amount of particles, the classical thermodynamic treatment is valid when the dimension of 
the confinement is much larger than the characteristic length of the system. Otherwise, the 
boundary of the system will cause obvious effects on the interior behaviors. 

For a binary fluid mixture, an important characteristic length scale is the correlation 
length, which will become infinite when the system approaches the critical state. Hence, the 
phase behavior of a confined binary fluid mixture can be affected much by the boundary 
of the system. Two common kinds of confinement are film and porous media, which cor- 
respond to the flat boundary case and the curved boundary case, respectively. The former 
is always modeled as a system confined between two slabs. For this case, when the system 
is quenched into the unstable region, the composition fluctuations with characteristic wave 
vectors parallel to the slabs can grow unrestrictedly while the ones with perpendicular wave 
vectors will be modulated by the confinement. As a result, the microphase structure in the 
system must coarsen laterally after the intermediate stage of the phase separation and the 
power law of the domain growth is different from the bulk case 1, 2, ^, 4, 5, 6|. In addition, 
each slab can be neutral or preferential to a certain component of the mixture so that the 
phase behavior can be coupled with the wetting phenomenon , . Compare to 



the flat boundary case, the porous media is much more complicated. The simplest model 
for this case is a binary system confined inside a cylindrical pore. Because there is only one 
free dimension for the cylindrical pore, the time required by the system to reach equilibrium 
is always so long that some metastable states are long-lived. The wetting behavior also has 

tinnnnnn 

great importance to the phase behavior for this case [12, ll3|, uM, Il5|, Il6|, ll7|, [ISj. 

Another interesting case is considered in this work. That is a polymer solution confined 
between two concentric cylindrical walls. This kind of confinement can be regarded as a 
combination of the two that mentioned in the last paragraph and hence is more general. 
When the radius of the inner cylinder converges to zero, the system will reduce to the one 
in a cylindric pore. On the other hand, if the radii of the two cylinders converge to infinite 
with a constant width between the two walls, the system will reduce to the one confined 
between two slabs. Thus, many results of this work can be easily generalized to the two 
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limiting cases. In addition, a polymeric system is always a good object for studying phase 
behaviors theoretically because the polymer chain length is a natural long enough length 
scale which makes the coarse-grained treatment and mean-field theory strictly valid. 

In the present work, we give a theoretical study on the phase behaviors of a polymer 
solution confined between two coaxial cylindrical walls. The self-consistent mean-field theory 
(SCMFT) is used to investigate the equilibrium structure of the system. The spinodal point 
of the system is determined by using the Gaussian fluctuation theory. The kinetics of the 
system are also studied by combining SCMFT with the Cahn-Hilliard theory. 



II. THEORETICAL FRAMEWORK 



Consider a polymer solution confined between two infinit long coaxial cylindrical walls, 
as illustrated in Fig. [T] The inner cylinder has the radius r and the outer one has radius R. 
The average volume fraction of the polymer is 0p, and that of the solvent is 0g = 1 — 0p. 
The canonical ensemble is used here, and the free energy of the system can be derived by 
the self-consistent mean-field theory 
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dr(/)p(r)</)s(r) + f drt/(r)(/)p(r) - ^ f dra;,(r)</),(r) 
T/0P^^.opiV_eQp_^- (1) 



In Eq. ([1]), Pq is the monomer density, which is defined as monomers per unit volume; 
a = P, S; X is the Flory-Huggins parameter, which quantifies the local interaction between 
each pair of polymer segments and solvent molecules; (/>a(r) is the volume fraction of species a 
at point r ; LJa{r) is the auxiliary field conjuncted to 0Q,(r); is the degree of polymerization; 

is the partition function of component a due to the kinetic energy, which can be regarded 
as a constant; is the partition function of a single molecule of component a. U{r) is 
the effective adsorbing potential to the polymers, which is assumed to be a rectangular- 
form potential for simplicity and is just put on the inner cylindrical wall in the present 
calculations. It is convenient to use cylindrical coordinates here and after. Then, U (r) has 
the form 

U{p,z)={ ° ^- (2) 
r + d < p < R 
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where Uq is a constant which quantifies the strength of the potential; d is the force range of 
the potential. To obtain the free energy, we need to solve a set of SCMFT equations, 



0P(P, Z) 
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dtq{p, z,t)q{p, z,N -t) 



i^piP, z) = X0s(p, z,t) + U{p, z) + ?7(p, z) 
^s(P^ ^) = X</'p(p, z, t) + ?7(p, z) 
(j)p{p, z,t) + 0s(p, z,t) = 1 



(3) 

(4) 

(5) 
(6) 
(7) 



Because the system has a rotational symmetry about the z-axis, the position-dependent 
functions in the SCMFT equations are independent on the polar angle. In Eq. ([3]), t is the 
coordinate along a polymer chain; the end-integrated propagators, g(p, z,t), satisfies the 
modified diffusion equation 



dqjp, z,t) ^ 
dt 6 
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q{p, z,t) - uJp, z)q{p, z,t) 



(8) 



with the initial condition, q{p, z, 0) = 1. All lengths in the present calculations are scaled 
by the Kuhn length b of the polymer, and thus 6 = 1 in Eq. ([8]). It should be pointed out 
that, to derive the above equations, the well-known Wiener measure is used to depict the 
conformation of each polymer chain. Hence our calculation should be restricted in the weak 
confinement case, i.e. R — r > Rg where Rg = (^NY^'^b. 

According to the Gaussian fluctuation theory [21], when^ there exists some fluctuations 
around the mean-field state, i.e. (/'o(r) = + (50a (r) ^], the free energy of the system 

can be written as an expansion 

F = + F(i) + + . . . (9) 

where F^^^ is the mean-field free energy illustrated by Eq. (P; F^^) = since the mean-field 
solution satisfies the SCMFT equations; F^^-* has the form as 



F(2) = ^ln<!det 
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In Eq. (fTOl) . the first term has no direct relations with current discussions; the functional 
^-'^^•'[{50}] in the second term is given by 



drdr'(C 



RPA\-1, 



r,r')50(r)50(r') 



(11) 



where S(P{t) = S<f)p{r) — 50s (r). The inverse of the RPA two-point correlation function in 
Eq. ([IID, (C^P^)-i(r,r'), is defined as 

(C^P^)-^(r, r') = C'\r, r') - xS{r - r') (12) 

where 

C{r, r') = C(r, r') - J dridr2A(r, r^)C-\v^, v,)A{r2, r') (13) 

In these expressions, the inverse operaters are defined through the relation 
/ driO~^(r, ri)0(ri , r') = 6{r — r'); the formulas of C(r,r') and A(r,r') can be found in 



ref. 



Many thermodynamic information of the system can be obtained from (C^^'^)~^(r, r'). 
Particularly, the condition that the smallest eigenvalue of (C^^^)"^(r, r') determines the 
spinodal point. In principle, the eigenvalues of (C^^'^)^^(r, r') can be derived by solving the 
eigenvalue problem 

I dr'(C^P^)-i(r, v')^nir') = A„vl/„(r) (14) 
where n denotes a complete set of quantum numbers. 

III. RESULTS AND DISCUSSIONS 
A. The case of neutral cy Under ic walls 

In this section, we mainly focus on the case of neutral cylindrical walls, i.e. Uq = 0. 
Because we always have R — r ^ b in this work, the boundary layers between the bulk of 
the polymer solution and the two cylindrical walls are less important. Thus it is reason- 
able to omit the boundary layers by letting q{p, z, t) satisfy the free boundary condition, 
i.e. = 0. After these treatments, it can be seen that the SCMFT equation set 

^ boundry 

has a solution for a homogeneous phase, in which (l)p{p, z) = 0p and (l)s{p, z) = 0g. This 
trivial solution represents the equilibrium structure of the system when x is small. How- 
ever, fluctuations will ultimately destroy the homogeneity of the system with the increasing 
of X- As mentioned above, the critical point at which the system becomes unstable (i.e. 
the spinodal point) can be derived from (C^^'^)~^(r, r'). Rather than solve the eigenvalue 
problem in Eq. (fT4|) directly, it is more convenient to expand (C^^'^)~^(r, r') with the eigen- 
functions of the Laplacian operator, V^. For the geometry of the system considered here. 
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(17) 



the eigenf unctions of are given by 
with the eigenvalues 

Ann.fc=(eH)' + fc2 (16) 

where Anmk is the normahzation factor; \m\ = 0,1,2,...; n = 1,2,...; A; is the one- 
dimensional continuous wave vector along z-direction. In Eq. (|T6l) . elT' can be zero or 
non-zero; when it is non-zero, according to the boundary condition, it is the n-th root of the 
following equation 

' J'a^l^'r) iV|'„|(eL-lr) 

where J\m\{x) and N\m\{x) are the |m|-th order Bessel function and the |m|-th order Neumann 
function, respectively. The radial part of the eigenfunction in Eq. ( |T5i) . Rnm{p), is 

{Sm,o ) eL ' = 

, f \m\ . JU^t^r) |„| |„| (18) 

|m| ^^'^ '' 

Expanding (C^^^)~^(r, r') by {ipnmk} leads to 

(C^P^)-i(r, r') = J2 anmkrnmk{r')^nmk{r) (19) 

nmk 

The expansion coefficient is given by 

1 . ^nmk fnn\ 
dnmk = + — X (20) 

where g{x) = ^— • As can be seen obviously from Eq. f|T9l) . the eigenfunctions of 
(C^^'^)^^(r, r') are also {ipnmk} with the eigenvalues ttnmk- Hence, the spinodal point can 
be derived from Eq. (l20l) that 

f 1 , ^nrnk 

Xsp = mm <^ + 



20g 4(Pp[l - g{XnrnkN)] 



--fi + ^l (21) 
Obviously, this result is the same as the non-confined case and is consistent with other 



related works 



22|. Physically speaking, the z-direction is a free dimension, thus the system 
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can undergo phase separations along this direction with no confinement, which is refflected 
by the confinement-independent spinodal point. 

It is worth noting that the mode corresponding to the minimum eigenvalue of 
r') is a constant, which corresponds a global translation of the system. This 
means that if x just equals Xgp, the system is in a critical state but with no phase separa- 
tion. However, once x exceeds Xsp; many fluctuation modes will be excited to induce the 
spinodal decomposition in the system. To study these modes in the early-stage of the phase 
separation, we firstly rewrite the free energy in Eq. ([1]) to a functional of 0p(r) by using the 
slow gradient expansion [23|, 

1 



dr 



X0p(l - 0p) + ^ In 0p + (1 - 0p) ln(l - (/)p) + — — (V0p)^ + 
iV 6b(p-p 



Fn + / dr 



(XSP-X)(0P-0P)' + T^(V0P)^ 



+ const 
(22) 



le other hand, the time 



where Fq is the free energy of the homogeneous state; the constant is omitted for simplicity, 
and the high-order terms are also omitted in the last step. On t 
evolution of (/>p(r) can be depicted by the Cahn-Hilliard theory 

5F 
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27| 



dr 



Hi 



(23) 



where r is a scaled time variable. By inserting Eq. (!22|) into Eq. (!23|) . it can be derived that 

1 



dr 



2(Xsp - x)H-i 



(24) 



where the fluctuation dcp-p = (f>p — (pp is used instead of (pp. 

Then we expand Scpp by {ipnmk}- Because of the axial symmetry of the system, it is only 
necessary to consider the fluctuation modes which are independent on the polar angle 6. 
Hence, 

5(Pp = J2 Cnkir)AnokRno{p)e''" (25) 

nk 

By inserting Eq. fl25!) into Eq. fl2^ . the evolving equation of the expansion coefficient, Cnk{T), 
is derived 



dc. 



nk 



dr 



'2(Xsp ~ X)^nOkCnk 



1 



\ 2 

''^nOk^nk 



The solution of Eq. fl26l) is 



c„fc(r) = c„fc(0)e""'=" 



(26) 



(27) 
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where 

180P 

= 2(X - Xsp) + k']-T^ Kr + k'T (28) 

When the system is thermodynamically unstable, x > Xsp- Then as can be seen from 
Eq. ( l28l) . for (e°)^ + A;^ < 360p(x — Xsp)' > 0, and the corresponding fluctuations 
increase with time. In particular, if 

(elf + e = 18MX - Xs.) (29) 

ujnk is maximized, and 

(^^nfc)max = 180p(x " Xsp)^ (30) 

This value corresponds to a maximum in the rate of increase of concentration fluctuations 
in the system. 

From the discussion above, it can be seen that the possible fluctuation patterns in the 
unstable system are dependent on 0p, r, R and A^. Generally speaking, these patterns have 
two parts, i.e. a radial part and a part along the ^-direction, as illustrated by Eq. ( l25l) . 
However, it is interesting that (ej)^ can be larger than 360p(x — Xsp) for many combinations 
of these parameters. For example, when (pp = 0.01, r = 3, /? = 50 and = 600, (e^)^ 
is not smaller than 360p(x — Xsp) until x exceeds 0.8. Note that Rq{p) = 1, then in the 
above cases, the fluctuation modes that can increase with time must be plane waves along 
the z-direction. This will lead to a periodic-like concentration profile along the 2:-direction 
in the early stage of the phase separation. 

Although the analysis above is in principle just applicable to the early stage of the phase 
separation in the system, it can indeed be concluded that the phase separation must be 
along the 2;-direction even in the late stage. This is because the system needs to put the 
interface vertical to ^-direction in order to minimize the interface area. As a result, even 
though the fluctuation mode may have a radial part in the early stage, only the ^-direction 
plane wave part can be maintained as time goes by. 
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B. The confined polymer solution with the presence of an adsorbing potential 



In this section, Uq in Eq. ([2]) is set to be positive. Hence, there is an adsorbing potential 
to the polymers at the inner cylindric wall. With the presence of this adsorbing potential, 
the boundary layer at the inner cylindric wall cannot be omitted and it needs to use the first 
kind boundary condition for q{p, z,t), i.e. q{p, z,t)\p=r = 0, because of the impenetrability 
of the inner wall. By contrast, the free boundary condition can still be used at the outer wall, 
0. Then the equilibrium structure of the system can be obtained by solving 
the SCMFT equations numerical 
direction implicit (ADI) method i 



p=R 

y. Especially, Eq. ([8]) is solved by using the alternating 



2S 



29|. 



In the present calculations, x is always set to be larger than Xsp ■ As such, one can expect 
some interplays between the adsorption of the polymer by the inner cylindric wall and the 
phase separation in the system. The main results are illustrated in Fig. [2l Note that the 
periodic boundary condition is used along the z-direction for convenience in the calculations, 
but the period along the z-direction in Fig. [2] is of no physical meanings. The strength of 
the adsorbing potential increases gradually from Fig. [2](a) to Fig. [2]^c). As can be seen from 
the figure, when the adsorbing potential is weak, it has few effects on the phase behaviors 
of the confined solution, and the system separates into a polymer rich phase and a solvent 
rich phase with the interface being vertical to the z-direction. However, accompanying with 
the increase of the adsorbing strength, a "wetting" layer can be formed by the polymer 
at the inner cylindric wall, as illustrated in Fig. Mjo). Furthermore, after the strength of 
the adsorbing potential exceeds a certain threshold, the interface between the polymer rich 
phase and the solvent rich phase becomes along the ^-direction and then the concentration 
profile of the system is independent on z (Fig. E^c)). 

As been discussed in Sec. IIIIA . the interface between the polymer rich phase and the 
solvent rich phase prefers to be vertical to 2;-direction in order to minimize the interface 
area. However, the ^-independent adsorbing potential, U, dislikes the inhomogeneity along 
the ^-direction. The competition between these two factors results in the variations in 
the concentration profiles, which are illustrated in Fig. [21 This phenomenon is similar to 
the "plug-tube" transition of a binary liquid mixture in a cylindrical pore. It is worth 
pointing out that the threshold of Uq, above which the transition happens, is always not 
very large, although it is dependent on the parameters of the system. In our calculation. 
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the value of this threshold is of order O-l^k^T per monomer. Hence if there exist certain 
strong interaction, such as the Coulombic interaction, between the cylindrical wall and the 
polymer, it is reasonable to omit the possible inhomogeneity along the ^-direction and the 
calculation can be reduced to one dimension {20 ]. 

Before ending this section, we would like to give some discussions on the periodic-like 
concentration patterns that are displayed in parts a and b in Fig. [2l Although the periods 
along the z-direction in parts a and b in Fig. [2] are of no physical meanings, these periodic- 
like concentration patterns can indeed be long-lived due to certain kinetic factors, which is 
observed by many experiments [3, 13] and simulations jlS, 14, 15]. It has been demonstrated 



in Sec. IIIIR that a periodic-like concentration profile along the z-direction can form in the 
early stage of the phase separation in the present system. In principle, this initial state will 
evolve to the equilibrium state by the process of domain coarsening. However, once the radial 
size of each domain reaches the size of the confinement, the coarsening process can only be 
achieved through the diffusion and the coalescence between adjacent domains along the z- 
direction, which will leads to a very slow kinetics. As a result, the periodic-like concentration 
pattern is long-lived, though is not the equilibrium state. One interesting point is that, if 
there is another faster process along with the domain coarsening, such as crystallization, then 
the periodic-like concentration pattern can be maintained. In reality, the inner cylinder in 
our model can be certain kind of cylindric adsorber and the outer cylindrical wall can be 
the interface between the adsorbed polymer solution and the outer environment. Hence the 
results presented here could be a hint to interpret the mechanism of the formation of the 
shish-kebab structure observed in the field of polymer crystallization [30, l31 ]. 



IV. SUMMARY 



In this paper, we give a theoretical study on the phase behaviors of a polymer solution 
confined between two coaxial cylindrical walls. The spinodal point of this system, which 
is derived by using the Gaussian fluctuation theory, is independent on the confinement. 
This is because of the existence of the free dimension, i.e. the 2;-direction. However, the 
fluctuation modes in the system is greatly affected by the confinement. Due to the kinetic 
analysis, the fluctuation modes in the early stage of the phase separation always have a 
component of a plane wave along the 2;-direction, which will lead to the formation of a 
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periodic-like concentration pattern along the axial direction of the system. The equilibrium 
structures are obtained by solving the SCMFT equations numerically and the interplay 
between the "wetting" phenomenon and the phase separation in the system is also observed. 

In particular, our results could give some hints to interpret the mechanism of the formation 
of the shish-kebab structure observed in the field of polymer crystallization. 
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Figure caption 

Figure [1] An illustration of the system considered. 

Figure [2] The concentration profiles obtained by solving the SCMFT equations with 
^ = 0.8, = 600, r = 3,R= 100, 0p = 0.1 and (a) Uq = 0; {h)Uo = 0.1; {c)Uo = 0.12. 
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FIG. 1: 
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